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I. INTRODUCTION 



In the present paper the possibility is explored that large-scale magnetic structures created during the 
QCD (Quantum Chromodynamic) chiral phase transition might lead to observable CMBR (cosmic microwave 
background radiation) correlations . This work is motivated by improved CMBR observations in progress, 
which promise polarization and temperature correlation measurements with multipoles I in the thousands. 
These magnetic structures could also be primoidal seeds of galactic and extra-galactic magnetic fields. A 
long-standing problem of astrophysics is the origin of the large-scale galactic and extra-galactic magnetic 
fields which have been observed. Of particular importance for cosmology is the possible seeding of these 
magnetic structures by primoidal, early-universe, magnetic structures. For a recent review see Rcf |l|. We do 
not, however, investigate galactic or extra-galactic magnetic structure seeded by the QCD phase transition 
in the present work, but center on CMBR polarization and gravitational wave correlations, 

In most of the theoretical treatments, including inflationary models, the magnetic fields arise from electri- 
cally charged particle motion. Considerations of nucleosynthesis and CMBR and the galactic magnetic fields 
have been used to constrain the magnitudes of the primoidal tangled (random) magnetic fields to values of 
about 10 -9 Gauss, although in a recent model [0 fields of two orders of magnitude larger are consistent with 
CMBR observations. Constraints on homogeneous primoidal magnetic fields from CMBR have also been 
determined ||. 

Early universe phase transitions are of great interest. The electroweak and QCD chiral phase transitions 
are of particular interest as common ground for particle physics and astrophysics. For the electroweak phase 
transition, using an Abelian Higgs model with the QED lagrangian included, magnetic field generation has 
been calculated ||[| in bubble collisions. For the QCD phase transition primoidal magnetic fields generated 
from charged currents at the bubble surfaces during the nucleation have been estimated |?],|| to be as large 
as 10 8 Gauss, and still be consistent with observed values of galactic and extra-galactic fields. There has also 
been a recent calculation Q of the CMBR power spectrum from density perturbations caused by promoidal 
magnetic fields that is similar to the calculation of Ref g), but for scalar perturbations. 

The present paper is motivated by our QCD instanton model of bubble walls formed during the QCD 
chiral phase transition |ifj|] and by the recent work of Forbes and Zhitnitsky, who have used a QCD domain 
wall |llj as a mechanism for generating magnetic fields which could evolve to large-scale galactic fields 
|fl2|| . Considering classical theory of bubble collisions, in which walls with the same surface tension as the 
colliding bubbles are formed within the merged bubbles, it was conjectured [|l0| that hubble-scale instanton 
walls might be formed, with lifetimes sufficient to form magnetic walls. In section II the description of the 
QCD phase transition in our instanton model and the possible resulting magnetic wall are discussed. It is 
shown that with an interior instanton wall, which is similar to QCD domain walls, the magnetic wall formed 
in the hadronic phase is similar to that modelled in Ref [03. In section III the evolution of the magnetic 
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structure to the time of recombination, and the CMBR polarization correlations are derived. In section IV 
the correlations from metric perturbations are derived, and in section V we give our conclusions. 



II. QCD CHIRAL PHASE TRANSITION AND MAGNETIC WALL 



In our model |10| of the bubbles of the hadronic-phase universe nucleating within the quark-gluon phase 
universe during the QCD chiral phase transition, which we assume is first order, we use the purely gluonic 
QCD Lagrangian density, C glue = G^G^/A. = d^A v - d v A^ - ig[A^A v ] is the color field tensor, 
defined in terms of the gluonic color field A^ — A"A n /2, where A™ are the eight SU(3) Gell-Mann matrices, 
with the notation that (p, v)= (1—4) and (a,b,...) = (1,2,3) are Dirac indices and color indices, respectively. 
Recognizing that it has been shown that QCD instantons can represent the midrange nonperturbative aspects 
of QCD, we use the instanton model (l3) for the color field 

p (x z + p ) 



r)- n Ap 2 

^liU W — / 2 _|_ ' ^' 

for the instanton and a similar expression with -n for the anti-instanton, where p is the instanton size and 
the 77" are defined as |1| 

Vij — e aij (3) 

with (i,j) = (1,2,3) and a,p,,v as defined above. The QCD instanton model, reviewed in Ref H|, has 
been a very useful representation of nonperturbative QCD for hadronic properties. For example, one can 
successfully predict the known scalar glueball candidates using this formalism 

Noting that the wall of the bubble separating the quark-gluon phase from the hadronic phase is gluonic in 
nature, and that the instanton model is successful in representating midrange nonperturative QCD, we have 
recently formulated |lC| ] an instanton model of the bubble wall. The starting point is the energy-momentum 
tensor of pure gluonic QCD: 

= G^Gl a - ± g ^G a /G a0a , (4) 

which gives in the instanton model at finite temperature in Minkowski space the spatial energy momentum 
tensor 

Tij ,inst = ( ViV ■ \ (5) 



(x 2 +p 2 ) 



and the energy density, 
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where p is the instanton size and N is the instanton density at the bubble surface at temperature T=T C , 
the temperature of the chiral phase transition, approximately 150 MeV. N is determined from the tunneling 
amplitude [|l3|,^4), and p has been found in finite temperature calculations [ fulfill to be p — 0.25/m. In Ref 
JTo| the energy density, T 00 was shown to be consistent with numerical calculations of the surface tension. 
T %3 can be used to calculate bubble collisions. 



2 



Recently, effective field models have been used to calculate QCD domain walls which could form within 
bubbles during the QCD chiral phase transition These domain walls have a space-time structure 

very similar to a wall composed of instantons, with the form given in Eq.(|l]). Such a wall could interact 
with nucleons in the hadronic phase to produce electromagnetic structures via effective magnetic and electric 
dipole moments [[l8| of the nuclcon field. This model was used jl2| to investigate possible primoidal magnetic 
fields at the time of the chiral phase transition that could lead to large-scale galactic magnetic fields. 

As pointed out in Rcf. JlQ[] , the collision of nucleating bubbles during the phase transition could lead to 
an interior gluonic wall. If the theory of classical bubble collisions can applied, the interior wall would be 
similar to the bubble walls, with the same surface tension. In other words, there would be an instanton wall 
with an energy-momentum tensor given by Eq.(||,|]). In a 1+1 model of colliding bubble walls based on QCD 
we recently found |l9| that, with instanton-likc boundary conditions, an instanton-like internal wall does 
seem to form in the collision region. Although this model is too simple for studying the nucleation problem, 
there have been a number of investigations of the QCD chiral phase transition using classical nucleation 
theories with effective Lagrangians |2^j2l|,|2^,^] . These studies find that in contrast to the electroweak 
phase transition where many bubles nucleate, collide and merge, the QCD chiral phase transition seems 
to proceed via inhomogeneous nucleation, with larger distance between bubbles and rather few nucleating 
bubbles involved in the transition to the hadronic phase. Therefore, it is likely that very few large-scale 
instanton walls were formed, and that the calculation of CMBR correlations from one magnetic wall, which 
is the picture used in the present study, is a good starting point. 

Recognizing that the mathematical form of the instanton and domain walls are very similar, the arguments 
of Ref. p2), including estimates of the lifetime of the interior QCD instanton wall, can be applied to estimate 
the primary electromagnetic wall that might have been formed at t~ 10 -4 sec. The magnetic wall is formed 
by the interaction of the nucleons with the gluonic wall, with the electromagnetic interaction Lagrangian 

where ^ is the nucleon field operator. This leads to the electromagnetic interaction with the nucleons 
magnetic dipole moment given in terms of the electromagnetic field tensor, F^ v by 

vint = ^rr*^75*F^, (8) 

and a similar interaction (without the 75) for the electric dipole moment, due to cp violation. From Eq.(^) 
one can estimate the magnetic field in the wall. In analogy to the classical theory in which the magnetic field 
within a magnetized object is given by the magnetic dipole moment density p4| ] within a factor depending 
on the shape of the object, the B-field is given by the matrix element of ^cr^j^^/ . For the gluonic instanton 
wall oriented in the x-y direction one obtains for B z = Bw = F 21 within the wall of thickness p 

A suppression factor of (pAqcd) 1 has been used in Eq.(^) since alligned dipoles tend to cancel, as discussed 
at length in Ref [fl2|| . The matrix element in Eq.(^) is estimated using the Fermi momentum in the plane of 
the wall during the QCD phase transition as Aqcd, giving < 1 I'(T2i75 v I / >= A qcd> with a factor of four 
from the spin and isospin degeneracy. Using Aqcd = 150 Mev, the resulting magnitude of the magnetic 
field at the wall is 

Bw - J^ A QCD, (10) 



which is essentially the same as the estimate of Rcf 12 , within the errors of the scale factors. The calculation 



of the electric field is similar, giving E z ~ B z ~ Bw — 10 17 Gauss (within the wall). 

Therefore our picture is that at the end of the QCD chiral phase transition there is a magnetic wall in the 
hadronic phase, which we model as 

B w (x) = Bffe-^^'e-^, (11) 
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or in momentum space 



where & _1 is of the scale of the horizon size, dn, at the end of the chiral phase transition (t ~ 10~ 4 s), 
b~ x = dn — few km, while M~ x ~ 0.2/m. Therefore, although By/ is very large, since the wall occupies a 
very small volume of the universe, such a structure is compatible with nucleosynthesis, galaxy structure and 
the present CMBR observations. 

In the present work we conjecture that the QCD bubble collisions lead to a magnetic wall given by Eq.(p"l|) 
at 10 -4 s and study the effects on CMBR polarization correlations and metric fluctuations. 



III. CMBR POLARIZATION CORRELATIONS 



In this section we investigate the polarization correlations arising from the electric and magnetic fields 
given in Eq.(^TJ). The primoidal magnetic wall of the present work is quite different from the tangled fields 
considered previously 26 
common with these studies. In addition to the polarization anisitropy which results from the magnetic wall 
itself, which gives rise to B-type polarization anisitropies discussed in the next subsection, in the subsequent 
subsection we also show that E-type anisotropies arise from scattering of the background radiation from the 
nucleon moments at the time of the phase transition, and show that the resulting power spectrum, CF E , is 
too small to be detected. 



, however, the evolution to the last scattering surface has many features in 



A. Polarization Correlations From Magnetic Wall 



In treating the temperature and polarization matrix for the Stokes parameters p7| , p8[ we use the angular 
representation of Ref. to derive the power spectrum of B-type polarization anisitropies, Cf B , which 

arise to a good approximation from the polarization source given by 

EP = J dnE(h,x, V )(Yi - V5Y] 1 ), (13) 

where the electromagnetic wave is propagating with k = fen, rj is conformal time, and the Y, m are the 

standard spherical harmonics. From our model the source at the wall is Eip — 5Bw/l2V%b 2 M n = Bw at 
the wall. The solutions of the Bolzman equation give the quantities Bf 1 fl29| at the time of recombination 
needed for the C BB . To get the power spectrum we must evaluate < B z (k, ?7)i? z (k', 77) >. Using the fact 
that exp(—k 2 d 2 H ) ~ 1.0 at the time of the chiral phase transition, 

< B z (k, v )B z (k', v ) > ~ B 2 w S(k x - k' x )5(k y - k' y ) < e -^/^ e -fcl 2 /4M= > 

^B 2 w d H e-W iM *5(k-k'). (14) 
This gives for the polarization power spectrum 

^.Sili^/tf^!, (15 , 

where the conformal time integral over the visibility function has been carried out and A77 is the conformal 
time width at the last scattering. The integral over the spherical Bessel function is carried out by using 
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the fact that ji(x) peaks at I and that the integral J dzjf(z) = 7r/(4Z) for large I. Therefore in the range 
100 < I < 2000 we have the approximate result 



BB 



1152M2At? 3 



(16) 



Using the parameters M n Ar] = 1.5 x 10 39 (from Refs. HQ), d H = 0.37 x 10 24 GW~\ and B w = 1.0 x 10 
Gauss, 
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Cf B ~ 4.25 x 10~ 8 Z 2 



(17) 



The result for the B-type power spectrum is shown in Fig. 1 by the solid line. There have been a 
number of investigations of the polarization predicted by inflationary models [j3l| which show that the B- 
type polarization in inflationary models is smaller than the E-type, and that 1(1 + 1)C BB peaks at Z-values 
about 100. Ref |5^] reviews, with extensive references, predictions of inflationary models and effects resulting 
from cosmological phase transitions. In Fig. 1 results from Ref |53| are shown. The curve shown by small 
circles is a string model normalized at I — 100 to our magnetic wall value, and with the same normalization 
the dashed curve gives a typical inflationary model result. Even the string model of Ref. [[33) is two to three 
orders of magnitude smaller than the T-spectrum at the peak, One can see that for I ~ 1000 the values of 
Cf B predicted by our magnetic wall picture exceed those of inflationary and topological models, and the I 
dependence is quite different. 



ln[l(l+l) cf B ] 



lOO lOOO 2000 / 

Figl. B— type power spectrum in magnetic wall model (solid), 
string (circles), inflationary (dashed) 



B. Polarization Correlations From Magnetic Dipole Scattering 



In this subsection we discuss another possible source of polarization correlation arising from the horizon-size 
instanton wall created during the chiral phase transition and leading to an orientation of nucleon magnetic 
dipole moments. As the nucleon moments are alligned along the wall the background radiation will scatter 
from these moments producing scattered waves and magnetic structure which would give polarization at the 
last scattering surface. Taking the thickness of the instanton wall as the scale for the volume associated with 
the magnetic dipole moments, The density of moments, n, is obtained from the strength of the magnetic 
field: 



B z = 



(18) 



w 
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where \i n is the neutron magnetic dipole moment and r\y is the wall thickness. Making use of the torque of 
the magnetic moment in a magnetic field B, dm./ dt = fi^s x B, with s the spin of the neutron, one obtains 
the magnitude of the scattered B-field at time t=10~ 4 s: B scat = n 2 ii 2 /2cos(8 s )kE exp(ikr)/r, with S the 
angle between the neutron spin and the direction of the incident radiation and Eq the E-field of the incident 
radiation. From this one obtains the Stokes parameter U, and in the notation of Ref p9] the polarization 
source 

B° 2 (k,rj) = ^(nul) 2 . (19) 

Note that E^fa-if) — B^ik^f), and results in E-type polarization anisotropies, C EE . 

Without a detailed calculation one sees that the resulting polarization correlations will be very small, since 
the extra k 2 dependence from the scattering from the moments introduces a factor of (l 2 /Ar/) 2 . Therefore 
the C EE polarization anisotropies resulting from the scattering from the nucleon moments at the time of 
the phase transition are too small to be measured. 



IV. METRIC PERTURBATIONS FROM MAGNETIC WALL AND CMBR POWER SPECTRUM 



In this section we derive the power spectrum from the gravitational waves arising from the magnetic wall 
of Sec. 2. The calculation is very similar to that of Ref. M, in which the power spectra was derived for tangled 
magnetic fields such as those considered in Refs. p5| , f26f , with various scenarios for the scale dependence. 
Although our model of the narrow hubble-size magnetic structure is quite different, we can make use of much 
of the formalism of Ref. || . 

The stress-energy tensor p4| ] for our model magnetic wall, which has only the 33 spacial component, is 

T 3 3 (k) = ^ J d 3 qB 3 (q)B 3 (k - q) (20) 

27T 3 0F , 



M, " w 



R p 8M £ 

D\yL 



at the time of the chiral phase transition. Note that we are assuming that the wall is of hubble size in the x-y 
directions and have omitted the parts of the expression for the k x , k y dependence shown in Eq. (|l^) . From 
Eq.(ECf) we obtain the B-wall power spectrum 



< T 33 (k,r,)T 33 (k',r)) > = j d 3 qd 3 q' < B 3 (q)B 3 (k - q)B 3 (q')B 3 (q)B 3 (k' - q') > . (21) 
In evaluating Eq.(|2l|) we use 



<e ^ S(k x )6(k' x )6(k y )d(k' y ) = e "^d^k - k')). (22) 

Using the notation of Ref Q with a tensor projection, 
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<T 33 (k,r ? )T33(k',7 ? ) > =^f(k 2 )d(k-k')), (23) 



with 



/ 2 ( fc2 ) = %d»Bw- (24) 

n 

With hij the tensorial perturbations of the Friedman universe, ds 2 = a 2 [—dr) 2 + (Sij + 2hij)dx l dx^]. The 
Einstein equations, using the representation hij = 2HQ^) of Ref ^9|, with V 2 Q-^ = —k 2 Q^ are 
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H + 2-H + k 2 H = 87rG*4/ 2 (fc 2 )- 
a a° 

The power spectrum for the tensor metric fluctuations are given by 

< hijP^hiiM > = 4|tf(k))| 2 <5(k-k')). 

Since Eq(p5[) was investigated in Ref 0, except with a different magnetic stress tensor, we use the solutions 
that they obtained. Assuming that the magnetic wall is formed at t=10~ 4 s in the radiation dominated 
epoch, at redshift Zi n , and neglecting perturbations created after the time of matter-radiation equilibration 
rj eq , for 77 > r} eq the approximate solution for H is 



(25) 



(26) 



H{k,rfi = A-KGif z eq ln{^) 



(27) 



Carrying out integrals over Bessel functions, the solution for the metric fluctuation power spectrum is for 
I » 1 



C, 



14 z . 

z, 



25~~ eq 



I'Vo I dz-J 2 {z/7 lo )J l+Z {z). 



Using the form of f(z) given in Eq(|24|), the integral in Eq(|2"8|) is approximately 

1 



dz ii! 2 { z / T io)Ji+z{z) 



8n 9 d H 4 0.106 



(28) 



(29) 



for M n r\ » I » 1. Taking t in at 10~ 4 s or z in = 1.17 x 10 s z eq and d H = 0.371 x 10 24 Gey _1 , the power 
spectrum is (I >> 1) 



l(l + l)Q ~ 6.9 x 1Q- 7 1 3 



(30) 



for B\y = 10 17 Gauss. Therefore, the metric perturbations from the QCD-induced wall result in CMBR 
effects large compared to other tensor perturbations that have been estimated for for I values of the order 
of 1000. 



V. CONCLUSIONS 



We conlcude that if the QCD chiral phase transition produces a gluonic wall of domain size and nucleon 
thickness at the time of the final collision of nucleating bubbles, and produces a magnetic wall as that found in 
the domain wall model of Refs 0,^2) , it would result in polarization correlations which have a /-dependence 
for I ~ 1000 different than other cosmological predictions and which should be measurable with the next 
generation of CMBR measurements. Also, there would be distinguishable temperature correlations arising 
from metric fluctuations. In order to investigate the collisions during nucleation to determine the details of 
the interior wall it is necessary to include quark and hadronic degrees of freedom to obtain the difference in 
the free energy in the two phases, a subject for future research. 
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